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We present an analytic approach to the estimation of beam asymmetry effects in CMB temperature 
and Unear polarization anisotropy experiments. We derive via perturbative expansions simple and 
accurate results for the case of an elliptical Gaussian window. Our results are applied to investigate 
the effect of beam ellipticity in the estimation of full-sky polarization correlation functions and the 
covariance matrix of power spectra. The relevance of this effect is also discussed by forecasting 
errors including beam asymmetry for current and future CMB experiments. 
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I. INTRODUCTION 

As high-resolution CMB experiments explore smaller 
fluctuations in the temperature anisotropy with high sen- 
sitivity, a better understanding of systematic effects is re- 
quired to make more accurate measurements. These sys- 
tematics have a direct impact on the abihty we have to 
improve the process of CMB data analysis at the level of 
map-making, power spectrum estimation and ultimately 
in constraining cosmological parameters. 

A common simplifying assumption in CMB data anal- 
ysis is to take the experimental beam response, i.e, the 
iso-contours of constant beam response, to have a per- 
fectly axisymmetric or circular shape with a Gaussian 
profile. This theoretical approximation introduces sys- 
tematic errors in the statistical analysis at angular scales 
comparable to the beam-width, a. Consistently, it bias 
estimates probing multipole orders ^ ~ 1 /cr in the spheri- 
cal harmonic analysis (i.e, the generalization of flat-space 
Fourier analysis for full-sky signals) of CMB experiments. 

As far as the main lobe is concerned, experimental 
beam responses for off-axis detectors are well-known to 
exhibit asymmetric shapes very well described by an el- 
liptical shape with a Gaussian proflle, as discussed for 
several experiments in the literature, e.g, Planck [||, p5[ , 
Maxima-1 Q and Python-V Q. However, the effect of 
beam asymmetry has been investigated only recently and 
the approach taken up to now has relied on semi-analytic 
p4[ or full numerical integration [||, pT[] . 

In this paper we shall introduce an analytic approach 
to address the problem of beam asymmetry in CMB ex- 
periments. In particular, we conveniently describe an 
elliptical Gaussian window in terms of a perturbative ex- 
pansion around a circular Gaussian one. As it will be 
shown below, this description allows a simple and in- 
tuitive discussion of the beam harmonic transform, the 
full-sky correlation and covariance matrices for both total 
intensity and linear polarization anisotropy observations. 

The paper is organized as follows: in §0 we present 
our analytic approach and derive the spherical harmonic 
transform of the total intensity beam. A detailed discus- 
sion of the effect of ellipticity to flrst order is provided 



[II. These results are validated numerically in §IV 



Results for linear polarization experiments are given in 
fjV|. We implement this formalism to calculate full-sky 



polarization correlation functions in §VI. Errors in tem- 
perature and polarization power spectra are discussed in 
Wn. Finally, we present a general discussion and our 
main conclusions in 



VIII. 



II. BEAM SPHERICAL HARMONIC 
TRANSFORM: TOTAL INTENSITY 

Let us consider the beam response, B, to the to- 
tal intensity sky distribution in a CMB temperature 
anisotropy experiment. For single-dish experiments with 
high spatial resolution, the beam geometry can be accu- 
rately described in the flat-sky approximation. Within 
this approximation, an elliptical Gaussian window func- 
tion can be expressed in cartesian coordinates. 



B{x,y) = Bo{aa,<Jb) exp 
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where we define aa and ah as the beam-widths in the 
major (x) and minor (y) axis, and the normalization is 
given by Bo{cra,(7b) = l/{2TT(Ja<Jb)- 

The Fourier transform of the flat-sky elliptical window 
is simply given by. 



B{kx,ky) = exp 



fc,v: 
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being kx and ky the modes along the major and minor 
axis of the ellipse, respectively. However, the Fourier 
analysis is only accurate for small patches of the sky {i.e. 
patches covering an area of a few deg^ or smaller) . 

For full-sky CMB analysis we shall introduce a decom- 
position of the window function in the spherical harmonic 
basis Yejn{0,4>), 
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where I w 7r/0 is the multipole order and him are the 
coefBcients of the harmonic transform, 
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being dfl = sm6d9d(p the differential sohd angle. Above, 
we have rewritten the elliptical window function B{9, (f>) 
in the (planar) polar coordinates, x = 9 cos(0 — uj) and 
y = 6'sin(0-cj). 
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where f{(p) = 1 — xcos^(0 — ui) describes the deviations 
from a circular (or axisymmetric) Gaussian window and 
the ellipticity parameter x = 1 — (ab/cTa)'^, is defined 
within the range f > x ^ 0- We have introduced an 
arbitrary phase uj which defines the orientation of the 
major axis of the elliptical beam in polar coordinates. 
The circular Gaussian window is thus recovered for the 
limiting case x = 0. 

However, the above integral for the spherical harmonic 
transform of the elliptical Gaussian window Eq.(Q) has 
no exact analytic solution and one has to resort to semi- 
analytic approaches or full numerical integration to eval- 
uate it (see eg, Q). 

In this section we shall show that Eq. {m can be solved 
analytically by introducing a convenient Taylor expan- 
sion of the elliptical (non-axisymmetric) window around 
a circular (axisymmetric) one. This perturbative expan- 
sion yields a series in powers of the ellipticity parameter 



B( 



= Bq exp 



Bo 



B{9)+x-^B(9) cos" {c^-Lo) 



+ 0{x^) (6) 



where the first term corresponds to a circular Gaussian 
beam B{9) = exp —9'^/2a^ of beam-width ab (the mi- 
nor axis of the ellipse; we shall denote ab = cr in what 
follows for simplicity) and Bq = [J dilB{9,4')]^^ is the 
beam normalization. 

The above expansion of the window function in real 
space Eq.(|y) leads to an analog expansion in harmonic 
space. 
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The n-th order term of the harmonic transform can be ex- 
actly integrated. In particular, only even m-modes have 
a non-vanishing harmonic transform n which reads (see 
Appendix 1 for the key steps of the derivation). 



*This is due to the azimuthal symmetry of the elliptical ; 
ometry, what is realized in the cos'^ (f factor of f{4>) in Eq.fl 
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being z ~ £^(T'^/2, Nim is the normalization of the spher- 
ical harmonics (see Appendix 1) and Lj^ (z) denotes 
the i'^^ order Laguerre polynomials of parameter a (see 



Eq.(AlS) for explicit forms). 

Replacing Eq.(|h into Eq.(Q) one gets the final expres- 
sion for the harmonic transform of the elliptical beam 



bem = cj"- A^/,„ e 
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where we define N^^ — N^-m/Bo, Bq = i?o/(27rcr^), and 
7^,,™ = {2iy + to)!/(22''+3™/2(;^ + m/2)\{i^ + m)\). Note 
that the circular Gaussian beam is recovered when x — 0; 
in which case, only the m — contribution is non-zero, 
bim = Nio exp[— ^^0-^/2] (5m, o- Combining the conju- 
gation rule for spherical harmonics, bg^ — (— l)'"6£-m 
(where non-zero m-contributors are even for an elliptical 
beam) and the reality condition of the beam transform, 
^irn — bim One sces that that both negative and positive 
modes have the same harmonic transform,, bi-.„i — bim- 
Therefore, in what follows, we shall assume m > with- 
out loss of generality. 

Eq.(0) is one of the main results of this paper. This 
equation demonstrates that the leading order correction 
to the circular Gaussian window from a given m-mode 
(for m > 0) is of order ©(x™^^). In other words, con- 
tributions from higher m-modes to the elliptical window 
function can be identified as higher order corrections to 
the circular Gaussian window. 

For high resolution experiments, cr <C 1 rad, the ellipti- 
cal beam harmonic transform is dominated by the axisym- 
metric or circular contribution to the window function, 
i.e. the to = mode. It is important to realize that the 
circular mode does no longer have a Gaussian profile, due 
to ellipticity corrections (see v ^Q terms in Eq.(||)). To 
leading order in the small x-expansion (j/ = in Eq.(ph), 
contributions from ttt, > modes are highly suppressed. 
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Therefore non-circular (higher-?7i) modes only have a 
non-negligible contribution to the harmonic transform 
with respect to the circular [m = 0) mode when a^^^ w 
1/x, which is well beyond the peak of the window func- 
tion. The peak of the window is determined from the 
leading order contribution to biQ Eq. (p^) . In fact, the 
peak location and width can also be accurately estimated 
from an effective circular Gaussian window of width, 
(7e//=a(l + x/4). 



^Higher-order corrections in the perturbative expansion 
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since x < 1 . Note also that the window function peaks at 
increasingly higher ^-multipoles as one considers higher 
(non-circular) TTi-modes. This is in agreement with recent 
numerical results p4| and provides a simple demonstra- 
tion for them. 



III. FIRST ORDER ANALYSIS: SLIGHTLY 
ELLIPTICAL BEAMS 

For most current and future experiments, such as 
Boomerang §, MAXIMA- 1 |6| and Planck |2|, the 
beam is only slightly elliptical, ie, the widths of the major 
(aa) and minor axis (a) of the beam differ by less than 
20%, 1.2 > (Ta/cr > 1 (x < 0.3). 

In this limit, a first order ellipticity correction to the 
circular Gaussian beam would give an accurate approx- 
imation to the actual beam harmonic transform, which 
yields for the modes m ~ and m = 2, 
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From this equation it is straightforward to see that for 
X ^ 1 the leading order contribution from the m ~ 2 
mode is a few percent of that from to = at the peak of 
the window a^ll^^^ ~ 1/2 (see Eq.(|l|)). 

However, notice that for the circular mode of the win- 
dow Eq. (jl^) the linear correction to the circular Gaussian 
window (second term in Eq.(n3)) is of the same order 
and peaks at the same multipole than the leading term 
in the non-circular (to — 2) mode (see Eq.(p^). There- 
fore, both corrections have to be included to compute the 
harmonic transform of the elliptical beam consistently. 
This is illustrated in Fig [^ for an elliptical beam with 
X = 0.3(cra/cr = 1.2) and ^fwhm = 10'. 



Eq.(y) ie, higher-order terms in v and m only modify sig- 
nificantly this definition for very large eUipticities x — 1- 
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FIG. 1. Ellipticity corrections to the harmonic transform 
for a slightly elliptical beam, with ellipticity x ~ 0.3 and 
resolution 6Ifwhm = 10'. (Solid lines) Leading order terms 
(for m — and m = 2). The first correction to the circular 
(m = 0) mode (dashed line) is of the same order and peaks 
at the same multipole than the leading order term for the 
non-circular (m — 2) mode (dot-dashed). 

Similarly, for highly elliptical beams, higher order x 
corrections to the circular mode become non-negligible 
and consequently higher non-circular modes have to be 
incorporated to calculate the harmonic transform accu- 
rately. Explicit expresssions for the window function up 
to second ord er in the ellipticity are given in the Ap- 
pendix A, Eq(A2C). This result arises naturally in the 
perturbative approach to the harmonic analysis of ellip- 
tical beams. 



IV. NUMERICAL INVESTIGATIONS 

In this section we shall validate the analytic results 
presented in the previous sections regarding the total in- 
tensity window. Although the same validation has been 
carried out for the linear polarization (see SQ), we shall 
concentrate here on the total intensity window as the re- 
sults for polarization are a straightforward generalization 
of the total intensity ones. 

First of all, we shall test whether the perturbative se- 
ries Eq.(||) is accurate and how fast it conver ges to the 
numerical solution. This analysis is done in §[VA. In 
addition, we shall see in §IVB to what extent the scaling 
solution for the higher ?7i- modes of the window, Eq.(|lO[), 
is a good approximation to the exact solution. As will 
be discussed in §IVC, prescriptions for an accurate com- 
putation of the window function from the perturbative 
solutions naturally lead to a criteria to how many higher 



(non-circular) ?Ti-modes have to be included in a consis- 
tent analysis of elliptical windows. 



B. Higher m modes contribution 



A. Probing the convergence 
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FIG. 2. Probing the convergence. Considering two differ- 
ent ellipticities of a given beam of width ^fwhm ~ W , we 
plot in the upper panel the n-th order expansion of bim and 
in the lower panel, the ratio of the n-th order correction to the 
0-th one. Both plots illustrate the nice and fast convergence 
of the X expansion. 

To test our approach, we shall compare the analytic 
results, Eq.(||), to a full numerical integration of Eq.(|j) 
using a Runge-Kutta method of fifth order |30| . The fast 
convergence of the analytical expansion is illustrated in 
Fig 1^. Indeed, in this figure we consider for one single 
beam size, ^fwhm = V81n2 a = 10', two different values 
of the ellipticity parameter, x '■ X = 0-'''5 {(Ja/cr = 2.0) 
and X — 0.30 {(Ta/a = 1.2) shown in the left and right 
panels, respectively. The upper panels display the ex- 
pansion fc^Q for various n. In both situations, the con- 
vergence is seen to take place for rather small n. Com- 
paring this two columns we see also that, as expected, 
the greater the beam ellipticity x is, the higher the num- 
ber of terms needed to reach the convergence. Both this 
statements will be precised quantitatively below. 

The lower panel illustrates this statement by drawing 
instead the ratio of the individual n-th order terms to 
the 0-th order one (for the same values of n than in the 
upper panels of Fig 0). 
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FIG. 3. Higher-7n modes in the elliptical window spherical 
transform. 

An analogous behavior as the one illustrated above for 
the m = mode convergence is seen for to 7^ 0. We 
remind here that odd m-modes are null and that only 
TO > modes are considered since negative modes have 
exactly the same harmonic transform, bim = bi-m (see 
^|)). Assuming this convergence, we now examine the 
amplitude of the higher m modes contribution, as they 
were derived analytically in §0. This is illustrated in 
fig.(0) for a beam of same width, i.e. 0fwhm = 10', and 
for 3 different values of x, namely x = 0-17 (ca/c — 1-1), 
X = 0.30 {aa/a = 1.2) and x = 0.41 {aa/cr = 1.3). We 
plot here the n-th order expansion of bim for m = 0,2,4 
where n is high enough so that this expansion is fully 
converged. 

The scalings demonstrated in Eq.(|l^) are clearly visi- 
ble. First we check that to > modes amplitude scales as 
(erf)™, making them not only sub-dominant (just a few 
percent contribution to the beam transform as compared 
to circular mode to = 0) but also shifting their peak to 
higher £ as to increases. Second we check also the scal- 
ing with the ellipticity, x™^^, that clearly implies that 
the smaller x, the more drastically the to > modes are 
suppressed. 

A direct comparison with the approximate scaling so- 
lution for the higher rn- modes Eq.(n9) is shown in Fig 
0. In particular, the plot shows both the fully converged 
expansion of these m. = 0, 2, 4 modes and the scaling so- 
lution for rather small values of the ellipticity x- We 
see that both the peak position and the amplitude are 
pretty well reproduced. Therefore the scaling solution, 



Eq(pi]|), is found to be a satisfactory description of such 
sub-dominant terms in the description of the eUiptical 
beam transform. Note however that the larger the beam 
elhpticity x, the worse this approximation turns out to 
be. 
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Table I. Required number of terms n in the ellipticity expansion, 
Eqteh, to achieve a precision greater than 1% up to £max = 5 ipeak 
for beams of different Sp-yvHM E'nd ellipticity x- 

The numbers presented in this table lead to another 
requirement. Indeed, as was discussed in § III the n-th 
order correction to the m = mode is of the same order 
than the leading order contribution to the m = n mode 
(to even) (see Eqs.(^ & (|lC|)). Thus to be self-consistent, 
the highest perturbative order, n, in the ellipticity cor- 
rections to the circular (m — 0) mode of the window, 
bgQ, should match the highest-m mode considered for an 
accurate computation of the full beam transform, bim- 
For example, the previous table implies that to handle 
properly the elliptical beam effects at a 1% precision till 
(^max = 5 (peak, e.g. for a beam of 6'fwhm = 10', we have 
to include m — 2 mode for x ~ 0-17 or 0.30 , while one 
has to include m — 2,4 modes for x — 0.40, and so forth. 



V. BEAM SPHERICAL HARMONIC 
TRANSFORM: LINEAR POLARIZATION 



FIG. 4. Probing an approximate computation of btrn for 
m = 2,4. Considering one beam of width Sfwhm = 10', we 
plot both the converged expansion of bf,„ and an approxi- 
mate evaluation of it as defined in Eq.([lO|). The left panel 
corresponds to x = 0.17 [pajo — 1.1), and the right one to 
X = 0.30 {jJa/u = 1.2). In both situations, the peak position 
and its amplitude are well reproduced. 



C. Prescriptions for an accurate analysis 

As was shown above, the convergence of the perturba- 
tive development is fast enough so that very few terms 
of the expansion are needed. To quantify this conver- 
gence and to define some useful prescriptions, we com- 
pare it to exact numerical integration and determine the 
order n of the biQ expansion needed in order to obtain an 
agreement better than 1% up to (.max — 5 £peak, where 
(peak denotes the maximum of the window function as 
defined i n Eq. (pll). Some prescriptions are summarized 
in table (IV C), where we write this order for two differ- 
ent beams of full width half maximum, ^fwhm = 5' and 
10' and a set of reasonable ellipticities. Even if we see 
that naturally the greater the ellipticity, the greater is 
the required n, as a matter of fact, in most of practical 
situations (see §VIII), 3 terms at most are needed. 

Note that this criteria is very stringent and if we re- 
quire, say, only a 2% accuracy at the peak level, only 1 el- 
lipticity correction is needed for x l£ 0-3, i.e. aa/cr < 1.2. 



The CMB radiation is expected to be linearly polarized 
as caused by Thomson scattering of CMB photons off hot 
electrons primarily at the surface of last scattering (see 
P2| , 0, |1^, [Q) while foreground Galactic emission is 
observed to be linearly polarized as well (see e.g. lO] 
for recent reviews and references therein). Thus we shall 
focus here on the detection of linearly polarized radiation 
and neglect circular polarization in what follows. 

The case for a linearly polarized beam with an elliptical 
shape can be treated in an analog way to the formalism 
developed in Sp^ for the total intensity beam. 

Linearly polarized radiation can be conveniently de- 
scribed in terms of the so-called Stokes parameters, Q 
and U (note that we use X to denote beam parameters, 
as opposed to sky parameters, X). Stokes parameters of 
a plane wave are related to the amplitudes of the electric 
field of the wave in two directions orthogonal to the wave 
propagation direction. 

Following standard notation (see e.g p3|, Q, p9J, 
pO{ , 1^ ) , the Stokes parameters of the beam are decom- 
posed in the spin-2 spherical harmonics basis ±2 y^m as. 
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Equivalently, the harmonic transform of a linearly polar- 
ized beam in terms of the so-called Gradient "G" and 
Curl "C" components reads, 
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from which it follows that 
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where the above expressions assume that the available 
power to each of the modes (G,C) is 1/2 of the total 
intensity (ie, we assume fully polarized detectors, with 
no sensitivity to circular polarization). Note that this G 
(Gradient) and C (Curl) components of the linear polar- 
ization are simply linked to the E and B ones, respec- 
tively, in the following way, 
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See ^^ for a pedagogical discussion of the E, B polar- 
ization modes. For a pure co-polar beam (ie, for an ideal 
optical system and telescope, see H), we have 



Q±iU = -B{ 



,±2i4 



(19) 



where B{9, (f>) is defined as in Eq.(H). Eq.(19) reflects the 
spin-2 nature of linear polarization in the (6, (j)) basis. 

Let us evaluate the harmonic transform of the linearly 
polarized beam. Using the parity symmetries for an el- 
liptical beam (see Appendix 2), 
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and the general (intrinsic to the definition of the G,C 
components) parity transformations, n. 



= h 
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P = G,C 



(21) 



one realizes that the harmonic transform of linear po- 
larization can be fully determined from one of the two 
components alone, say G. Moreover, both negative and 
positive modes have the same harmonic transform. Thus, 
in what follows, we shall assume rn > without loss of 
generality. 

In full analogy with the total intensity computation 
(see §0) we introduce a perturbative expansion of the 
elliptical beam. 
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This expansion can be exactly integrated for any order 
in an analogous way to the total intensity case and yields 
(see Appendix 2 for the basic steps of the computation) , 



■'■Note that these conjugation rules are consistent with |2d] 
& IM, and are inconsistent by a factor (—1)™' with respect 
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(23) 
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where we define N^^ = -f.^"" Mem/{'iV2Bo), and the 
coefficients ji,^m-2 are the same than those defined for 
the total intensity Eq.(||). The normalization Mim is 
given in Appendix 2 along with the basic notation for 
the spin-2 harmonics. 

Note that, up to the normalization N^^, the linear po- 
larization beam transform Eq. (kSiI is formally the same 
than the total intensity one EqT(@), with the index for 
the m- modes shifted by to — )■ ttt, — 2. This shift is in- 
troduced by the difference in spin-index s between the 
linearly polarized beam s = 2 (see Eq.(n9) above) and 
the total intensity beam s — 0. 

In particular, Eq.(p3) shows that the m = 2 mode dom- 
inates the harm,onic transform, of the linearly polarized 
elliptical beam. Note that for a circular Gaussian beam 
(x = 0) only the to = 2 mode is non- vanishing. To lead- 
ing order in the ellipticity expansion (1/ = in Eq.(23)) 
one finds that contributions from m > 2 modes are sub- 
dominant. 
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(24) 



with 777, = 4, 6, . . . In particular, one finds the same sup- 
pression of higher to- modes with respect to ?7i = 2, in 
full analogy to the results for the total intensity Eq. {nm . 
Also, the linear polarization window function peaks at 
increasingly higher £-multipoles as one considers higher 
TO,-modes, as was the case for the total intensity window. 
The expressions for the first non-zero TO-contributors 
(m = 2 and m — A modes) up to the first ellipticity 
correction are 
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and therefore. 
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(27) 



and the same expressions hold for negative modes, as 
bf_„^ — bf.^. This is in full analogy with the scaling 
relation between higher TTi-modes found for the total in- 
tensity window Eq.(O). 

Whenever the beam ellipticity is fairly large, one has 
to consider higher-order corrections in the ellipticity to 
compute accurately the window function. Explicit ex- 
presssions for the window function up to secon d ord er in 
the ellipticity are given in the Appendix A, Eq(B22). 



VI. FULL-SKY POLARIZATION CORRELATION 
MATRIX 



Qeff = 2^[Dl,M,0,u;)]* bf*, af„^ 



(31) 



Linearly polarized radiation is described by the total 
intensity T and the Stokes parameters Q and U. If the 
CMB polarized radiation is Gaussian distributed, one 
needs, a priori, 6 statistical quantities to characterize 
correlations among them. It is more convenient to use 
linear combinations of the Stokes parameters with differ- 
ent parity properties, the so-called E (or Gradient G ) 
and B (or Curl C) modes, for which only 4 correlations 
are non- vanishing. Namely the correlation between E 
and T modes and the 3 autocorrelations. 

Following [|42| we will consider the correlation matrix 
M between 2 arbitrary measurements in the sky 
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M(ai,S2) = 
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U1U2 
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(28) 



where 1,2 denote the directions ni,n2 in the sky. The 
cross terms < T1U2 >=< Q1U2 >= as required by 
symmetry under parity transformations (see e.g. |20[| ). 

The entries of the correlation matrix are defined as 
follows: 



<PlP2> = <P*ffinuLUi)Peff{n2,LU2)> 



Peffin,Uj) 



dn D{(l),9,u}) P* P 



(29) 



where Peff is the result of convolving the polarized beam 
P = T, Q, U with the sky, P = T, Q, U. 

In this formalism the "scanning strategy" of a given 
experiment is obtained by specifying the Euler angles 
as a function of time t, {(l){t),9{t),uj{t)), where n = 
n{9(t), 4>{t)) gives the pointing direction of the beam and 
U!{t) is the rotation angle around the pointing direction 
n which specifies the orientation of an asymmetric beam 
(eg, the major axis orientation for an elliptical Gaus- 
sian beam) with respect to a fix reference orientation (eg, 
w(i = 0)). 

Accordingly, the rotation operator D{(j),6,uj) acts on 
the beam so that it takes all possible orientations with 
respect to a fix reference frame in the sky |4y] , |y . Simple 
scanning strategies allow a convenient decomposition of 
the rotation matrix D{(p{t),9{t),oj{t)) for the implemen- 
tation of fast methods to compute the full-sky convolu- 
tion 1^21 7 o- In what follows we shall suppress the time 
dependence of the Euler angles to simplify notation. 

Decomposing the polarization field in spin harmonics 
one finds the following expressions for the Stokes param- 
eters of the beam convolved with the sky (see Eq.(5) of 
@ and Eq.(39) of |l). 
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(32) 



where we define, af^ = {a2jm + a-2,^m)/2\/2 and 
*aFm = (a2,fm - a-2,£m)/2%/2. 

Note that for the linear polarization parameters Q, U, 
an overall factor of 2 accounts for the fact that both G 
and C modes contribute equally to the transform of an 
elliptical beam. 

The polarization correlation matrix can be easily com- 
puted making use of symmetry properties of the rotation 
matrices and the addition theorem of rotations (see Ap- 
pendix 3), 



< T,T2 >= E brcf 



<TiQ2>-2E bJ^C, 



TG 



<TiU2 >cx C/'- =0 



< Q1U2 >oc Cf"- = 



<0iQ2>=4E bf'^C^i 



<UiU2 >-4E bf^Ci 



(33) 

(34) 

(35) 
(36) 

(37) 

(38) 



where we have used the fact that Cf'~"' — Cf-^ = as 
follows from the general property that the T, G harmonic 
coefficients of a field transform differently under parity 
than the C harmonics. The power spectra are defined 



< O^m (^e'm' >= CfSil'Smm' , P = T, G, C (39) 

and we have introduced the "2-point window functions" , 



,pp' _ v^ ^i ,p, ,p' 

MM' 



T, G (40) 



where hereafter we drop the tilde to denote beam quan- 
tities, {i.e. we take P ^ P) to ease notation. We remind 
that the first non- vanishing contributions to the total in- 
tensity bj^ = bim enter at tti = 0, while linear polariza- 
tion beams bf^ have the first non-zero contributor from 
m = 2. 



The rotation matrix D\,j,j^j above reads (see Eq.(2) in 
§4.7.1 and Eq.(5) in §4.7.2 of HI) 

(41) 

where the Euler angles (a, /3, 7) of the resuhing rotation 
matrix are (Eq.(6) in §4.7.2 of ||; see also g]), 

cot a = cos 02 cot((/>i — ^2) — cot^i sin 6*2 esc ((/)i — ^2) 

cos /3 = cos 61 cos 62 + sin 6*1 sin 6*2 csc((?!)i — 02 ) 

cot 7 = cos 01 cot((/)i — 02) — cot 02 sin 02 csc((/)i — 4)2) (42) 

and the orientation of the beam at pixels 1 and 2 is given 
by cji and 0^2 • The irreducible rotation matrices (see 
Eq.(2.17) of and Eq.(4) in §4.3.1 of |§), read 






,(/?) = }_^(-l) 
t 



{£ + ny.il - n)\{e + my.{e - m) 



1/2 



tl{£ + n- ty.{£ - m - t)!(t + m - ny. 
X (cos/3/2)2^+"-"-2*(sin/3/2)2*+'"-" (43) 

where t is summed up for all values which yield non- 
negative factorials. These matrices give the dependence 
of the polarization correlation functions on the separation 
(or lag) angle in the sky P = rii ■ n2. Thus, predictions 
for the polarization correlation matrix for cosmological 
models of the sky signal convolved with an elliptical win- 
dow assuming a particular scanning strategy are given by 
Eqs(|),(|2|),@ & (H). 

The polarization correlation matrix thus obtained can 
be used to compute the likelihood functions and the 
Fisher information matrix for a given sky realization of 
a cosmological model convolved with an elliptical beam. 



A. Slightly Elliptical Beams 

Provided the beam ellipticity is small (x <C 1) a first 
order ellipticity correction to the circular Gaussian beam 
yields a good approximation to the elliptical beam trans- 
form (sec §111). Consistently, one can expand the polar- 
ization correlation functions to first order in the elliptic- 
ity expansion. 

For this purpose, first we need to explicit the first terms 
of the 2-point window, 



h^^' - D^ h^ h^' 



n^ hP hP' 



(44) 



where we have taken the real part of the 2-point function 
as we want to compute the polarization correlation func- 
tions which are observable and therefore real. Eq.(|4^) 
for the case of the total intensity [P = T) is in agree- 
ment with Eq.(33) in ||34| . Introducing the expansion of 
the (1-point) windows &?* to first order in the elliptic- 
ity Eqs(|l2|),(|l|),(|5|) & (|6|) into Eq.(||) one obtains the 
following expressions for the correlation functions. 
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where the sum involving the rotation matrices D'j^j ^^, 
restricted to ^ > M + M', and 
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+ (-1)^44 (tt - /3) cos(2a - 47) + cos(4a - 27) (49) 



where we have to replace a — > a — W2 , 7 — * 7 — ^i to 
introduce the beam orientation in the above equations, 
as determined by Eq.(pl|). We have used the symme- 
try properties of the d-rotation matrices (see §4.4 Eq.(l) 

in 11) dMM'{l3) = i-l)^''-^dMM'm,dM-M'iP) = 

{—l)^'^d,MM'{'^ — /?), and the reahty condition on the 
correlations functions. The irreducible d<„„-matrices can 
be expressed in terms of Legendre Polynomials by relat- 
ing the previous to the spin-s harmonics (see Eqs(3.4), 
(3.11) in |l|), 



TT Correlation Function 
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For s = 0,2 one gets (see Eqs(A4) & (B2)), 
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and _2^™(cos/3) are given in Eq.(B4). This explic- 
itly shows that the ellipticity (asymmetry) of the win- 
dow function introduces a dependence of the correlation 
functions on the scanning strategy, as parameterized by 
the angles (a,/3,7). We stress that the above equations 
are appropriate for a full-sky analysis, since the small- 
angle approximation is only taken for the beam geometry 
(which is of small extent in radians). 

In Fig ra we display the temperature correlation func- 
tion < TT > (we drop sub-indices labeling sky pixels to 
ease notation) for a slightly elliptical beam Eq.(p3). We 
assume a flat power spectrum, Ce — Constant. This al- 
lows us to emphasize the effect of the window function ir- 
respective of the underlying cosmological model assumed. 
For the case shown in Fig a, we assume that the telescope 
scans the sky in ecliptic latitude {i.e. a = j — 0) and 
that the beam hits a given sky pixel always with the same 
orientation {i.e. we consider correlated pixels for beams 
aligned in the sky) which provides an upper limit to the 
effect of ellipticity on the correlation functions. This is 
because scanning strategies which observe a given sky 
pixel with a different orientation of the beam each time 
it scans over it, tend to average out the impact of beam 
asymmetry on full-sky estimators. 
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FIG. 5. Temperature correlation function for a slightly 
elliptical beam. It assumes a flat power spectrum, 
Ci — constant. (Short dashed line) Correlation function for 
a circular Gaussian beam. (Long dashed) Linear Ellipticity 
correction assuming the beam scans the sky in ecliptic lati- 
tude and a fixed beam orientation in the sky. (Solid) Total 
correlation function (Gaussian -I- linear correction). Thick 
lines assume x = 0.3{(7a/<y ~ 1.2) while thin lines correspond 

to X = 0.17((Ja/(J = 1.1). 
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FIG. 6. Same as Fig M but for the absolute value of the cor- 
relations in log scale to emphasize small residual correlations 
induced by the beam at large angular separations. Only lines 
for an ellipticity of x = 0.3 are shown for clarity. 
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FIG. 7. Cross-correlation Temperature-Polarization. Con- 
ventions are the same than in Fig |5[ 
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FIG. 9. Linear Polarization correlation function in terms 
of the Stokes Q parameter. Note that for the case shown (a 
flat power spectrum), < QQ >—< UU >. Conventions are 
the same than in Fig H. 
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As seen in the plots, the hnear elHpticity correction 
to the circular Gaussian window introduces an anti- 
correlation for pixels separated by < era (the major 
axis beam-width). This is because pixels within this an- 
gular separation are seen as a single smeared pixel. Al- 
ternatively, for 9 > CTa the ellipticity increases the cor- 
relation between sky pixels. This correlation peaks at 
9 = 2a a, where it yields a 20% correction (for x — 
0.3) to the Gaussian correlation function and decreases 
monotonously for larger distances, as expected. How- 
ever, a closer look (see Fig O) reveals that some small 
residual anti-correlations induced by the ellipticity (at 
the level of 10^^) remain at large distances which might 
be a reflexion of the limitations of a linear order analysis. 
Non-linear terms in the ellipticity expansion are expected 
to cancel out these long-range correlations. 

Cross-correlations for temperature-polarization and 
linear polarization auto-correlation functions are shown 
in Figs (0),(|i) & (p|) , dlOl) , respectively. In particular, we 
see that the effect of ellipticity is comparable for the tem- 
perature auto-correlation < TT > and cross-correlation 
functions < TQ > (at most a 15% correction to the Gaus- 
sian correlation for x = 0.3), while tends to be less sig- 
nificant for the linear polarization < QQ > (just a few 
percent correction). Notice that the angular scale for 
the transition between negative and positive ellipticity- 
induced correlations (first and second bumps in the log- 
scale plots) is shifted in the correlations involving linearly 
polarized windows with respect to the case for the tem- 
perature (or total intensity) discussed above (see Figs (g) 

&(!)). 



VII. COVARIANCE MATRIX 

In this section we shall discuss the covariance of the 
power spectra for elliptical beams in the presence of un- 
correlated|j noise, following the standard formalism de- 
veloped for circular windows |2|] (see also Q , ||2y] , Q , 
Q, Q & Q). In particular, we shaU use this formal- 
ism to estimate error-bars for the power-spectra for ellip- 
tical window functions. For this purpose we shall assume 
that the circular mode of an elliptical window, which has 
a non-Gaussian profile, yields an approximately unbiased 
estimate of the actual error-bars, as we shall argue below. 

The covariance of the temperature power spectrum, 
Cj can be easily computed for the circular mode [m = 0) 



of the window function [pT[ 



A(CJ) 



(2£ + l)Mfsky 
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(52) 



Similarly for linear polarization, one obtains from the 
lowest TO-mode contribution {m = ±2) ^^, ||44| , 
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where P — G, C, and for the cross-spectra 
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(54) 



where the factor of 2 in the polarization windows ac- 
counts for equal contributions from m = ±2. The fac- 
tors A£ and fsky in the above expressions account for 
the binning in £-space used (we assumed Ai = 75 for 
all experiments) and the fraction of the sky observed 
by the experiment, respectively. The weight per solid 
angle is w = (a'i^^u}pix)~^ while the noise per pixel 
c^pix = s/ ^tpix depends on the detector sensitivity s and 
the observing time per pixel tpi^. The pixel solid angle 
^pix — ^FWHM X ^FWHM- The above expressions for the 
noise associated to power spectra estimation, Eqs.(|52|)- 
(p4) assume that all detectors in the experiment have 
the same noise properties and main beam response. 

Note that polarization power spectra have twice as 
much noise per pixel as the temperature spectrum since 
only half of the total power is available to each polariza- 
tion mode (G and C). This is accounted for through the 
normalization of the window functions (see factor v^ in 
Eq.(|6l))J3 

We stress that the above expressions only include the 
leading order in the m-mode expansion of the elliptical 
window. However this is a good approximation to (i.e. 
, it is the dominant term in) the exact window function 



for elliptical beams as discussed in §111. In principle, this 
analysis could be rigorously extended to include higher 
771- modes of the window by computing the aim's of the 
sky map convolved with the elliptical window, from which 
the power spectra of the convolved map and their as- 
sociated errors can be calculated (see appendix A. 2 in 
||l6|). However thorough numerical analyses ( |^l|; see 
also S) show that an azimuthally symmetrized compo- 
nent of the window yields an unbiased estimate of the 
power spectrum within a few percent, what suggests that 
non-circular modes of the window function can be safely 
neglected, at least for slightly elliptical beams. 



Hhe noise is assumed to be uncorrelated between differ- 
ent pixels and between temperature and linear polarization 
measurements. 



**Alternatively, one can define different pixel weights w for 
temperature and linear polarization, w^ = 2 uF (see e.g. fcol) 
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Predictions for the theoretical error-bars for the tem- 
perature power spectrum for current CMB experiments 
are shown in Fig O. Experimental parameters are taken 
from |l|], H (MAXIMA-1), (H), |l[ (Boomerang), |] 
(Archeops) and [p6[ (Planck). The figures used corre- 
spond to averages among channels and they only intend 
to be illustrative. Note that for Archeops and Planck 
the experimental numbers given are just nominal. It is 
seen that the pixel noise blows up the error bars at mul- 
tipoles i ~ 1000, except for the Planck satellite experi- 
ment. Main differences between error forecasts for differ- 
ent experiments are due to the sky coverage and noise per 
pixel (for a single channel). It is also observed that the er- 
ror bars computed for a Gaussian window underestimate 
those of an elliptical window (computed according to 
Eq.(n2|). However to first order, the error bars for an el- 
liptical beam can be well approximated by using an effec- 
tive Gaussian window of the form &£o = exp[— £^(7aCr/2], 
where CTq and a arc the major and minor axis of the el- 
lipses of constant beam response. 



ERRORS IN POWER SPECTRUM (SINGLE BOLOMETER) 
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FIG. 11. Errors in the power spectrum estimation for cur- 
rent and future experiments. It assumes an underlying stan- 
dard ACDM model. Pairs of lines (above and under the mean 
power spectrum) show different estimates for the theoretical 
error bars according to different choices of the window func- 
tion: short-dashed lines are predictions for a circular Gaussian 
window of beam- width given by the minor axis of the ellipse 
o, solid lines correspond to the elliptical window to first order, 
Eq.(|l3), while long-dashed lines are obtained from an effec- 
tive circular Gaussian window of beam-width (Je// = ^OaO , 
being Oa the major axis of the ellipse. In all cases, a binning 
in £-space is used of width A£ = 75. 

A detailed analysis of the expected error bars in 
the power spectra estimation including polarization for 
Planck (single 100 GHz channel) is summarized in Fig n2. 



As discussed above, the high sensitivity of Planck allows a 
clean recovery of the CMB power spectra up to ^ > 1000 
with a single channel data (except for the C-polarization 
mode, see below). In fact, pixel noise starts blowing up 
the error bars for the temperature anisotropy power spec- 
trum at ^ > 1500 (see upper panel in Figfi2|). For the 
cross-spectrum temperature-polarization (G-mode) one 
finds that pixel noise becomes dominant at £ > 1000 (see 
middle panel) whereas for polarization (G-mode), this 
happens at lower multipoles !l « 1000 (see lower panel). 
We have checked (although is not shown in Fig |l2| for 
clarity sake) that for the polarization C-mode error bars 
become pixel-noise dominated at ^ < 500 as the signal is 
typically (i.e. for standard CDM models) found at a few 
percent level of that in the G-mode. 

Beyond these multipoles [i.e. for smaller scales) the 
effect of the ellipticity of the window becomes signifi- 
cant. Moreover, using a circular Gaussian window clearly 
underestimates error bars for elliptical beams approxi- 
mately computed according to Eqs.(|l2|) & ( p5| ) for the 
total intensity and linear polarization windows, respec- 
tively. 
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PLANCK 100 GHz (SINGLE BOLOMETER) 
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FIG. 12. Errors in the power spectrum estimation from a 
single 100 GHz detector of the Planck satellite. It assumes 
the same experimental parameters than those given in Fig [III 
(bottom right panel). Solid lines show the mean theoretical 
power pectra and their error bars for realizations of the sky 
convolved with a Gaussian beam. Dashed lines show the ana- 
log error bars for the case of an Elliptical beam. Upper panel 
displays the temperature anisotropy power spectrum, middle 
panel shows the cross temperature-polarization (in terms of 
the G-mode), while the bottom panel corresponds to the po- 
larization (G-mode). 
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VIII. DISCUSION 

As cosmic microwave background (CMB) experiments 
image the sky at finer spatial resolution with higher sen- 
sitivity, new relevant systematic effects have to be prop- 
erly taken care of in the process of data analysis in order 
to consistently extract cosmological information down to 
the smallest scales probed by the experiment. The asym- 
metry of the beam response is becoming an increasingly 
important issue which has been largely neglected until 
recently in CMB studies. 

In this paper we have introduced an analytic approach 
to describe the effect of beam ellipticity in CMB experi- 
ments. This approach is based on a perturbative expan- 
sion around the geometry of a circular Gaussian beam, 
which yields a series expansion of the elliptical Gasu- 
sian beam in powers of the ellipticity parameter. There 
are several advantages of introducing a perturbative ap- 
proach to discuss beam ellipticity: 

• It provides a simple and convenient way of inte- 
grating the beam harmonic transform for the total 
intensity and linear polarization. 

• In most of current experiments the beam ellipticity 
is small (we shall refer to these as "slightly elliptical 
beams" in what follows), i.e, the beam fuUwidths 
along the major and minor axis differ by 10-20 % 
at most. This implies that, in practice, the pertur- 
bative expansion truncated to low orders describes 
the harmonic transform with high accuracy up to 
very high multipoles. 

• The perturbative expansion allows a simple quali- 
tative discussion of the role that different m-modes 
play in the beam transform (see e.g, section III). 
In particular, the relative weight of these modes is 
assessed by working out how they depend on ex- 
perimental parameters (e.g, width and ellipticity 
of the beam). This information cannot be directly 
extracted from a non-perturbative solution. 

• The full-sky polarization correlation matrix can 
be most simply discussed for the "slighly ellipti- 
cal beams" for which deviations from the circular 
Gaussian beam results can be explicitly derived. 

In particular, we have obtained analytic solutions for 
both the total intensity (temperature anisotropy) and lin- 
ear polarization window functions. The main results are 
given in §|| Eq.|) & §0, Eq.(|2|). 

Our findings show that the circular (m = 0) mode 
dominates the total intensity window function, although 
the first non-circular (higher m) modes can not be ne- 
glected in a consistent analysis. The reason for the lat- 
ter is that higher m-modes in the beam transform can 
be identified amongst the higher-order corrections in the 
ellipticity expansion around the circular Gaussian win- 
dow. This provides a simple explanation for previous 
semi-analytic and numerical results in the literature. 



For linear polarization, we found that m — 2 is the 
dominant mode but again, higher modes (m = 2, 4, . . .) 
must be included to compute accurately the window func- 
tion. 

Numerical integration validates our approach and pro- 
vides practical prescriptions for how many terms in the 
perturbative expansion of the circular mode of the win- 
dow have to be taken to achieve a given accuracy. This 
in turn translates directly in how many non-circular 
(higher-TTT,) modes contribute non-negligibly to the win- 
dow function of the elliptical beam (see § [V A ) . 

We have implemented our analytic solutions for the el- 
liptical window function to derive expressions for the full- 
sky polarization correlation functions for elliptical beams 
(see §|Vl|). In particular, we have derived simple ana- 
lytic expressions for slightly elliptical beams, taking into 
account the beam orientation and scanning strategy of 
a given experiment. We find that, for simple scanning 
strategies, the ellipticity of the beam induces additional 
correlations of the order of 20% for small angular sep- 
arations (few beam-widths) with respect to a circularly 
symmetric Gaussian beam. 

Finally, we have investigated the impact of beam asym- 
metry in error estimation for CMB power spectra in the 
presence of uncorrelated noise. We find that error bars 
for a circular Gaussian window largely underestimate 
those of an elliptical window when the pixel noise be- 
comes dominant. However, a good approximation to the 
actual error bars is given by an effective circular Gaus- 
sian window of beam-width Ueff = \/o\fa, being aa and 
a the major and minor axis of the ellipse. Note that, 
for slightly elliptical beams (x ^ 0)i fe// « cr(l -t- x/4) 
which is approximately the width of the circular mode 
of the elliptical window, &£0i as discussed in §||. This 
explains why for quasi-circular windows, Ueff — y^o^, 
provides an accurate estimate of the power spectra error- 
bars. 

We shall emphasize that, in our approach, we introduce 
the experimental beam in the time stream^ while the " ef- 
fective beam" in the pixel domain is the result of multiple 
observations of the same sky-pixel with different orienta- 
tions of the beam (and possibly with a different noise 
level), for general scanning strategies. This implies that 
non-circular modes of the " effective beam" are expected 
to cancel out to some extent and therefore the " effective" 
circular component of the beam should yield an almost 
unbiased estimate of the C|, as shown by recent numer- 
ical analysis (see Wu et al 2001). Therefore the nomi- 
nal ellipticity in the time domain will be typically larger 
than the final "effective" ellipticity on the map. In the 
discussion of t he e stimated errors in the power-spectrum 
presented in §VII, we take the "effective" ellipticity to 
be the same than the one defined in the time stream and 
thus our estimates must be taken only as upper limits to 
the actual effect of window ellipticity. 

The issue of beam asymmetry here discussed is partic- 
ularly relevant for future high-resolution and sensitivity 
CMB anisotropy experiments, especially those measuring 
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also polarization, such as the Planck satellite. 

In a future work |12|, we shall validate the elliptical 
model for the beam asymmetry presented here in the 
presence of other systematic effects (non-elliptic beam 
distortion/asymmetry, pointing errors, other sources of 
noise, etc). Such analysis will show under which circum- 
stances beam ellipticity is a major systematic effect in 
a realistic analysis of a CMB experiment. Some recent 
work along this lines has already been done for the Planck 
satellite [|0|, although the formalism used is only valid 
for small-patches of the sky. 
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APPENDIX A: PERTURBATIVE EXPANSION 

OF THE ELLIPTICAL BEAM HARMONIC 

TRANSFORM 

In this appendix we present the key steps for the 
derivation of the spherical harmonic transform for the 
total Intensity beam Eq.(^). 

In the flat-sky limit {9 -^ 0) the elliptical beam shape 
can be expressed in cartesian (x,y), 



B{x,y) = Bo{aa,(Jb) exp 



X 



y 

2al 



(Al) 



where we define Ua and ct^ as the beam-widths in the 
major x and minor y axis, and the normalization is given 
by Bo(cra,crb) = l/(27raoCr{,). 

For analysis on the sphere, it is more convenient to in- 
troduce (planar) polar coordinates to describe the beam, 
X = 9 cos 4> and y = 9 sin 0, 



Bi 



Bq exp 



2al 



fW 



(A2) 



where /((/>) = 1 — xcos^ describes deviations from the 
circular (or axisymmetric) Gaussian window and the el- 
lipticity parameter x = 1 — (ab/iTa)^, is defined within 
the range 1 > X ^ 0- The circular Gaussian window is 
thus recovered for the limiting case x — 0- For the sake 
of simplicity, we have taken the beam to be pointing to 
the north pole of the sphere {9 = 0). 

The spherical harmonic transform of the total intensity 
beam is defined as, 



dn B( 






(A3) 



where, dfl = d9 sin 9d(f), and spherical harmonics are de- 
fined as, 



Yi,ni9,(t>) ^ NirnPricos9) e'"''^ 



Ni„i = 



2£-Ki {e~my. 
47r y {e + my. 



(A4) 



(A5) 



where P™ are the Legendre polynomials and the spher- 
ical harmonics obey the conjugation property Y^*„ = 
(-l)"y£_™. Replacing Eq.(R|) in Eq.(0) we get. 



bern — (— l)'"^£-m 

X d9 fim9Pf"'{cos9) d(PB{9,(j)) e'""*' (A6) 

^0 Jo 

In the flat-sky hmit (6* < 1 rad, ^ > 1), 

Pf"'{cos9)^r"'J^{i9) (A7) 

where Jm is the m-th order Bessel function of the first 



kind. In this limit, the above integral Eq.(A6) reads 



d9 9J,n{e.9) / d<f)B{9,(l))e'' 
Jo 



(A8) 

In order to solve this integral analytically, we introduce 
a convenient perturbative expansion of the beam in real 
space in powers of the ellipticity parameter x, 



B( 



r H 

Bo B{9) exp^x ^ ^os^ , 



*«(''i:(|?) 

ri=0 



„2n J, ^ 



(A9) 



being B{9) = exp —9^/2a^ | , which yields an analogous 
series in harmonic space. 



E^ 



(")X 



r,(0) __ 



+ fc^lix + 0(x') (AlO) 



The perturbative expansion Eg. ([A9|) allows to factor- 
ize the 2D-integrals of the beam harmonic transform 
Eq. ([A^) in two ID-integrals for 9 and </> respectively. 
Thus the n-th order term of the beam transform can be 
expressed as follows: 



uM 



= (-l)"W,_,„4:]ifW 



(All) 



with. 



4:^=r"(2a2)-" 



^(») 



d9 9 



2n+l 



M£9)B{9) (A12) 



^^^2n I —i2m(b 

COS (be ^ 



(A13) 
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Making use of Eq.(6. 631.1) of H and Eqs(13.1.27), 
(13.6.9) of one gets, 



r(") 



2+m 



{n~m/2)\ 



2™/2 



^n-m/2(^) 



with z = fa'^/2 and, 



^(«) 



27r 



2n! 



(n + m/2)!(n-m/2)! 



(A14) 



(A15) 



'(") 



for m even, and Xm = for m odd. The fact that odd 
TO-modes do not contribute to the harmonic transform 
is due to the parity symmetries of the eUipse. Thus the 
rt-th order term of the expansion Eq.(All) is given by, 



b^l = 



m 22n+m/2 



^2+m -z T(m) 

(n + m/2)! %>-W2V 



(A16) 



which replaced in Eq.( A17 ) yields the final expression. 



where Niq = yj2i + l/47r, and negative modes [i.e. m — 
—2, —4) have to be included as they have the same har- 
monic transform than positive modes, i.e. bi^m — bim- 
Similar expressions for the harmonic transform to first 
order in the ellipticty are given in §[11, Eqs.jl^) & dl3). 
Note that for a circular Gaussian window, x = Oi one 
gets as expected, beo — A'^foexp[— ^^cr^/2] and be„i = 
for m y^ 0. 



APPENDIX B: PERTURB ATIVE EXPANSION 

OF THE LINEARLY POLARIZED ELLIPTICAL 

BEAM HARMONIC TRANSFORM 

The aim of this appendix is to provide a detailed 
derivation of the harmonic transform for linearly polar- 
ized elliptical beams Eq.(p3|). The spherical harmonic 
transform of a linearly polarized beam can be written in 
terms of the Stokes parameters Q & [/, 



'Jim 



a N f 



oo 

e-^ ^7.,™4'"n2)x'^+™/' (A17) 

where we define iV/^ = Ni^m/BQ, Bq — Bo/{2'Ka^), and 
7^,„ = (2i/-hm)!/(22''+3™/2(j,_,_TO/2)!(iy + m)!). The first 
Laguerre Polynomials are. 



6« = — 

'" 2V2 

6^ =— ^ 

'™ 2V2 



dn 
dn 



[Q - iu) 2y;„, + {q + tu) -2y;, 
(Q - iu) 2r;„,, - (Q + lU) -2Y,\ 



(Bl) 



-M 



(m). 



L'P{z)^l ■ LP(z)-m + l-z 

4"^^) = i[(TO+l)(TO + 2) + z(-4-2m + z) 



where we define the spin-2 spherical harmonics as [j 



Lt\z) 



Yf.n,.=Mt 



±2J-lm 



Im ± 



2P"\cOSt 



^imrj) 



where 



(TO+l)(rn + 2)(m + 3) 

-3(to + 2)(to + 3) + z(9 + 3m - z) 



(A18) 



Me, 



' (/-2)! 
(/ + 2)! 



Ni, 



(B2) 



(B3) 



and higher-orders can be obtained from the recurrence 
relation (see Eq.(4.18.1) in (2|) 

1 



and we define a generalization of the Legendre polyno- 
mials for spin-2 harmonics, Fj 



Li-'Hz)- 



{2n -1-X + m)Li"!\(z) - (n - 1 + m)L^,^l\{z) 

(A19) 



±2 



pr(cos( 



£ — m^ 



i- 

V SI 



sin^ei 



e{e~i))Pl^{cost 



In most practical situations the beam ellipticity is rather 
small, X <C 1. In these cases, one only needs to compute 
the first terms (2 or 3 terms account for the beam trans- 
formup to v ery la rge multipoles with high accuracy, see 



+ (£ + to)^^P;!!i(cos( 



sm 



T 



IV q . Table |IVCD 



For example, the beam harmonic transform up to sec- 
ond order in the ellipticity expansion has non-vanishing 
contributions only from the modes m = 0,2 and 4, which 
read, 



('(^- 1) cos 6IP™ (cos 61) 
- {e + m)P^i{cos9)') (B4) 



sm 



The above quantities obey the following parity condi- 
tions: 



beo = Neo e 



-t'a'^ 12 



X^2^2 






-1 c 



(B5) 



011= N in — I (J e ' 



x(i-^^V 



X tf4_4 „-f^cr^/2 



bi^^NiQ ^ta" e 
128 



(A20) 



^^Note that ±2^^™ = W^^m ± *^<m, according to the nota- 
tion used by IM. 
**The ±2-P/' polynomials are simply related to the Gim in 

M, m, [g: ±2Pr(^) = G+„(a.)TG7„(x). 
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which imply that, 

±2^/m — 



(€ + m)! 
(€-m)! 



M/ 



^2^£m 



-2iin4> 



(B6) 



(B7) 



_2Fr(cos0) « -(-1)"£™+V„_2 



(B13) 



and thus, 






what allows us to recast Eq.(Bl) in a more convenient 
way, 

1 



4V2 



(B14) 



h^ - 



—i 



dn 
dfi 



iQ~lU)2Yern + {Q + iU)-2Yi, 

{Q~iU)2Yi„,-{Q + iU)-2Yt, 



-2imct> Introducing the ellipticity expansion Eq.(|A2|), one can 
solve the integral to any perturbative order. 



-227710 



uGir^ 



1 



(B8) 



2V2 



Af,7,7 12 K, 



(«) ^-(77) 



(B15) 



For a pure co-polar beam (ie, for an ideal optical system 
and telescope, see B), we have 



Q±iU = -B{e,, 



±2id 



(B9) 



where B{9,(j)) is defined in Eq.(A2). We have assumed 
that the beam response is measured in the co- and cross- 
polar basis defined on the sphere, ctco and Ccross, accord- 
ing to Ludwig's 3rd definition 



Ceo = Sm (p Gg + COS (j) (70 
Across = COS (j) (Tg- siu (j) U^ 



(BIO) 



where ag and a^ are the usual spherical polar basis. Such 
CO- and cross- polarization basis, Eq.(BlO), is obtained 
by parallel-transporting the local cartesian basis defined 
at the north pole, CTx and cry, along great circles through 
the poles of the sphere (see e.oj^ for a discussion). 

Replacing Eq.(B£) into Eq.(BS), one sees that the first 
term in Eq. ( |B8[ ) is non- vanishing only for negative m- 
modes while the second term is non-zero for positive m- 
modes alone. What is more, the parity properties of the 
G and C modes, 



be„, = I 



G 



P = 



n = -* bf_ 

G,C 



(Bll) 



imply that the harmonic transform of linear polarization 
can be fully determined from one of the two components 
alone, say G. Moreover, both negative and positive modes 
have the same harmonic transform. Thus wc shall assume 
TO > below with no loss of generality. In this case the 
harmonic transform of the G-mode is simply given by. 



h^ - 



M, 



£. m 



2V2 



dnB{9,(l))^2Pricos0) e 



7(777-2) 



(B12) 



In the flat-sky limit (6* < 1 rad, £ > 1) 



^^This corrects the expression for t he small-angle limit in 
[ p5[ : the pre-factor £™+-^ in Eq.(B13) corrects the pre-factor 
r^^™ inEq.(4.32) of S. 



with, 



fin) 



pn+2 



IiZ = ^^i'2^) / d9e'''+'J^.2i^9)Bi9) (B16) 



(.27r 

K^"'' = d(P cos^" (P e 

Jo 



7(777-2) 



(B17) 



Noting that the above integrals are basically t he sa me 
t han those for the total intensity beam Eqs(A12) & 
( |A13| ), but replacing m by to. — 2 everywhere, they can 
be integrated in the same way, 



-(„) _ (n- to/2 + 1)! (™_2) . X 

-'^777— " om/2 ^ -^77-777/2+1 V^/ 



2777/2 

with z = £2 0-2/2 and 

Kin) ^ 2^ 



2n! 



2277 (r7, + m/2-l)!(n- to/2 + 1)! 



(B18) 



(B19) 



for even modes ?n > 2, and K, 



(77) 



for TO, odd. There- 
fore the n-th order term Eq.(B15) in the beam expansion 
is given by. 



lG(77) 



27r 



:Mj 



2nl 



22n+m/2 " "\n + m/2 - 1)1 



-z r(™-2) 

^77-m/2+l 



a'" e ~ L'll" ^',^ , ^{z) 



(B20) 



which introduced in Eq.(B15) flnally gives. 



if _ ni-2 prG „-2\^^ r{m-2)(^ v+m/2-l 



u=0 



(B21) 



where we define iVf^ = -f"" Mim/{4:V^Bo), and the 
coefficients 7,^,777-2 are th e same than those defined for the 
total intensity Eq.(A17), except for the subindex which 
is 771 — 2 here instead of to there. 

For most of actual experimental beams the ellipticity 
is rather small, x <C 1. As discussed in Appendix A, 
a second order analysis of the beam ellipticity is already 
very accurate to very large multipoles as comparison with 
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numerical integration shows (see j |lVC , Table IVC for 
specific prescriptions depending on experimental param- 
eters). Thus, expanding the beam harmonic transform 
to second order in x one gets non- vanishing contributions 
only from the modes m — 2,4 and 6, 






1 - ^^^2 



+ 



X 



16 



Of 4 



2V2 



X ^2^2 ^-.VV2 



1 + X 1 



-£^a^ 



,G ^ _^ X^ 
" 2V2 128 



^4^4 ^-Pa-/2 



transform to first order in the ellipticty are given in §M, 
Eqs.(p5|) & (|26|). Note that for a circular Gaussian win- 
dow, X = 0, one gets, bf^^ = -(A^f o/2\/2) exp[-^2^V2] 
and bfjj^ = for \m\ > 2. N ote that, as argued above (see 
paragraph und er Eq .( B17 )), the linearly polarized beam 
transform, Eq.(B22), can be straightforwardly obtained 
from the total intensity beam transform, Eq.(A2C), by 
replacing in the latter m by m — 2, and including a mul- 
tiplicative normalizing factor of — 1/(2\/2) appropriate 
for linear polarization modes, see Eq.(|Bl[). 
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